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Abstract. The centrality dependence of dynamic fluctuations of the transverse 
momentum and the net charge can signal the approach to local thermal equilib- 
rium in nuclear collisions. I explore this signal by comparing transport-theory 
calculations to STAR and PHENIX data at a range of energies. In particu- 
lar, I find that this model can describe PHENIX data on the dependence of 

■ fluctuations on the p t range in which they are measured. 

o 

■ Keywords: Event-by-event Fluctuations, Relativistic Heavy Ion Collisions 
O ! PACS: 25.75.Ld, 24.60.Ky, 24.60.-k 



O 



- 1—1 

X 

5-H 



1. Introduction 



Measurements of the event- by- event fluctuations of the mean transverse momen- 
tum exhibit substantial dynamic contributions [[I]- PHENIX and STAR data for 
Au+Au collisions show that such fluctuations increase as centrality increases [E3 E]- 
In addition, the measured dynamic net-charge fluctuations may show an analogous 
behavior Importantly, the event-averaged transverse momentum (pt) - a quan- 
tity unaffected by fluctuations - exhibits a strikingly similar behavior [0J E). 

I ask whether the onset of local thermal equilibrium can explain the common 
centrality dependence of these quantities. I seek a nonequilibrium explanation be- 
cause the strongest observed changes occur in peripheral collisions, where equilib- 
rium is not expected. This is important because similar behavior can accompany 
the hadronization of quark-gluon plasma, but only if the matter is thermalizcd [ 6 . 
My discussion follows the transport-theory formulation in ref. [Cj. Here, I present 
work in progress confronting new measurements of the variation of pt fluctuations 
as a function of the observed p t region see fig. [21 

Dynamic fluctuations are generally determined from the measured fluctuations 
by subtracting the statistical value expected, e.g., in equilibrium [|Hj. For particles 
of momenta pi and P2, dynamic multiplicity fluctuations are characterized by the 
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robust variance 
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dp 1 dp 2 r(p 1 ,p 2 ), 



(1) 



where (• ■ ■) is the event average. This quantity depends only on the two-body 
correlation function r(pi, p 2 ) = iV(pi, p 2 ) — N(pi)N(p 2 ), see also [|JJ. It is obtained 
from the multiplicity variance by subtracting its Poisson value (N) . To make Raa 
robust, we divide by (iV) 2 to minimized the effect of experimental efficiency [0. 
For dynamic pt fluctuations one similarly finds 



where Spu — pu — (pt)', STAR measures this observable. The observable measured 
by PHENIX satisfies F Pt as N (SpaSptz) /2cr 2 when dynamic fluctuations are small 
compared to statistical fluctuations a 2 = (p 2 ) — (pt) 2 - 

Thermalization occurs as scattering drives the phase space distribution within 
a small fluid cell toward a Boltzmann distribution, which varies in space through 
the temperature T(x, t). The time scale for this process is the relaxation time v . 
In contrast, density differences between fluid cells must be dispersed by transport 
from cell to cell. The time needed for diffusion to disperse a dense fluid "clump" 
of size L ~ (|Vn|/n) _1 is ~ vL 2 /v 2 h , where v t h ~ 1 is the thermal speed of 
particles. This time can be much larger than v^ 1 for a sufficiently large clumps. 
Global equilibrium, in which the system is uniform, can be only obtained for t 3> td- 
However, the rapid expansion of the collision system prevents inhomogeneity from 
being dispersed prior to freeze out. 

Dynamic fluctuations depend on the number of independent particle "sources." 
This number changes as the system evolves. Initially, these sources are the indepen- 
dent strings formed as the nuclei collide. The system is highly correlated along the 
string, but is initially uncorrelated in the transverse plane. As local equilibration 
proceeds, the clumps become the sources. The fluctuations at this stage depend on 
number of clumps as determined by the clump size, i.e., the correlation length in 
the fluid. Dynamic fluctuations would eventually vanish if the system were to reach 
global equilibrium, where statistical fluctuations would be determined by the total 
number of particles. 

2. Mean p t and its Fluctuations 

In [[7] I use the Boltzmann transport equation to show that thermalization alters 
the average transverse momentum following 



where S = e~ is the probability that a particle escapes the collision volume 
without scattering. The initial value {pt)o is determined by the particle produc- 
tion mechanism. If the number of collisions J\f is small, 5* ~ 1 — M implies the 





(Pt) = (Pt)oS + (Pt)e(l - S), 



(3) 
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Fig. 1. Mean p t (left) and its fluctuations (right) vs. centrality with data [EJ El- 
Preliminary A data is from [0] • Centrality in all figures is determined by the number 
of participants relative to the b=0 value. 



random- walk-like increase of (p t ) relative to (pt) - For a longer-lived system, energy 
conservation limits (pt) to a local equilibrium value (pt)e fixed by the temperature. 

As centrality is increased, the system lifetime increases, eventually to a point 
where local equilibrium is reached. Correspondingly, the survival probability S in 
(PJ) decreases with increasing centrality. The average p t peaks for impact parameters 
near the point where equilibrium is established. The behavior in events beyond that 
centrality depends on how the subsequent hydrodynamic evolution changes (pt) e as 
the lifetime increases. Systems formed in the most central collisions experience a 
cooling that reduces J3J with proper time t as (pt)e oc r -7 [|7j. 

Thermalization can explain the behavior in fig. ^ The survival probability is 
S — exp{— J" F v(t)<1t} w (tq/tf)™, where v = (crw rc i)n(r) for a the scattering 
cross section, ?; r ei the relative velocity, and to ; _f the formation and freeze out times. 
Longitudinal expansion implies n(r) tx t , yielding the power law with a = isqTq. 
To fit the measured centrality dependence, I assume a = 4 and 7 = 0.15 in central 

1 /2 

collisions, and parameterize S(N pairt ) by taking a oc N pll] . t and rp — tq oc N paTt , 
where A^ par t is the number of participants. I take the same a for all species, as 
appropriate for parton scattering. 

Dynamic fluctuations depend on two-body correlations and, correspondingly, 
are quadratic in the survival probability. In [d I add Langevin noise terms to the 
Boltzmann equation to describe the fluctuations of the phase space distribution. 
A simple limit is obtained when the initial correlations are independent of those 
near local equilibrium, (5p t i5p t2 } = (SptiSp t 2) S 2 + (6pti5p t 2} e(l - S) 2 ; this form 
was used in [ [7] . Alternatively, if the initial correlations are not far from the local 
equilibrium value, I find 

(SpaSpa) = {SptiSp t 2)oS 2 + {6p tl 6 Pt2 ) e (l ~ S 2 ). (4) 

Here I use (@J, which provides somewhat better agreement with the latest STAR 
data in the peripheral region where thermalization is incomplete (and my model 
assumptions most applicable). As before, the initial quantity (5pti5pt2)o is deter- 
mined by the particle production mechanism, while (5pti5pt2) e describes the system 
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Fig. 2. Dynamic pt fluctuations with PHENIX data for central Au+Au as a func- 
tion of the acceptance 0.2 GeV < p t < pt, max [li- 



near local equilibrium. 

To estimate (6pti8pt2)o for nuclear collisions, I apply the wounded nucleon 
model to describe the soft production that dominates {pt) and (SptiSpt2), to find 

Wtxvnio N ^ \ 1 + Raa J W 

[d- The pre- factor is expected because measures relative fluctuations and, 
therefore, should scale as N~^. t . The term in parentheses accounts for the normal- 
ization of to (N(N — 1)); Raa scales as N~^. t ISR measurements imply 
(Spti5pt2) PP / (pt)p P ~ 0.015. I use HIJING to estimate R pp and Raa- 

Near local equilibrium, spatial correlations occur because the fluid is inhomo- 
geneous - it is more likely to find particles near a dense clump. These spatial 
correlations fully determine the momentum correlations since the distribution at 
each point is thermal. The mean p t at each point is proportional to the temper- 
ature T(x), so that (SptiSp t 2)e ~ / ^(xi,X2)<5T(xi)5T(x 2 ), where r(xi,x 2 ) is the 
spatial correlation function, ST = T — T, and T is a density- weighted average. I take 
n (oc T 3 ) and r to be Gaussian with the transverse widths Rt and £t, respectively 
the system radius and correlation length. In [|7] I obtain 

(S P nSpa)e = (6) 
1 + Raa 

where Raa is given by [ 7 . The dimensionless quantity F depends on £,t/Rt- I 

i /2 

compute F assuming R t oc N p ' alt -. F = 0.046 for £t/Rt = 1/6. I again use HIJING 
to estimate Raa- Calculations are compared to PHENIX data in fig. G] (right), and 
agree within the large uncertainties. 

PHENIX also reports F Pt in the acceptance region 0.2 GeV < p t < Pt, ma x 
as a function of pt,max, with the intent of identifying the role of jets in produc- 
ing fluctuations To exhibit the acceptance dependence of F Pt , I use JQ and 
(0) to estimate (5p t iSp t 2) ~ <t 2 Raa, where a 2 = (p 2 ) — (pt) 2 - This result is ex- 
act for short range transverse correlations. I then write the observed quantity 
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F Pu0 bs « ((N)(SptiSp t 2)/2a' 2 ) ohs oc (N) ohs R AA:0hs . The observed multiplicity 
grows strongly as pt jmax is increased, since (iV) bs is the integral of the pt dis- 
tribution over the acceptance. To determine the pt.m&x dependence of Raa. obs, let 
the probability that an individual particle falls in the acceptance region be a. The 
average number of detected particles is {N) b s — a(N). For a binomial distribution, 
one has (N 2 ) obs = a 2 (N 2 ) + a(l - a){N), so that (Ql implies R,AA,obs = Raa, the 
value for full acceptance. It follows that 

r pt, Obs ~ r p t \ I ) 

where F pt and (N) are for full acceptance. I emphasize that this estimate only 
includes trivial acceptance effects. I omit contributions to the pt, ma x dependence 
from the onset of jets, as well as changes in the correlations due to restricting p t . 

The estimated pt jma x dependence of F Pt accounts for the data in fig. [3 The 
solid curve is computed for (iV) bs oc 1 — e^ Pt ' xnax ' T (l +pt,max/T) for a thermal 
distribution at the freeze out temperature obtained from the measured (pt)\ the 
dashed curve for T = 600 MeV is included to show the T dependence. These 
curves are arbitrarily normalized, since my computed values in fig. ^ fall above the 
most central PHENIX data (but below STAR, see fig.[3J|. Jets may account for the 
discrepancy at the highest max- 

The agreement in figs. ^ and [5] is excellent - and therefore somewhat puzzling 
since I have neglected collective flow [ HOj . Flow must contribute to fluctuations at 
some level. However, thermalization and radial flow are manifestations of multiple 
scattering in different regimes; the flow velocity is uniquely defined only when local 
equilibrium is established. Both effects enhance pt, so that the flow contribution 
may be partly concealed by my parameter choices. 

Preliminary STARp t and net charge fluctuation data in fig.[3]show fluctuations 
at 20, 130 and 200 GeV beam energies. To compute this energy dependence, I 
assume a oc N, Raa oc N" 1 and take N from data. The pt fluctuation data show 
a tantalizing similarity to the calculations, although the experimental uncertainty 
must be reduced to firmly establish the behavior at high multiplicity. Net charge 
fluctuations in fig. El (right) are characterized by Vd yn = R ++ + R — 2R^ for 

R ab 

given by [ I12j . I therefore expect Vdyn to satisfy i^dyn — 
use pp data to fix v$ and compute v £ /vq assuming the longitudinal narrowing of the 
correlation length seen in balance function data [^2- The agreement of in-progress 
calculations with preliminary data [OH I12| is encouraging. 

Deviations of the calculated fluctuations from data in central collisions at the 
highest energies are expected due to flow [ and jets [ ^ . Such deviations may 
be seen in fig. [3] Nevertheless, I stress that hard scattering effects - including 
jet quenching - are incorporated in my calculations via the HIJING Raa in I©- 
HIJING alone does not describe the data discussed here [lllj. 

Thermalization is consistent with the acceptance dependence in fig. [21 I sur- 
mise that (JZJ holds for a range of models, since my argument was rather general. 
That said, precise measurements may yet reveal jet and flow contributions H0| - 
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Fig. 3. Dynamic pt (left) and net charge (right) fluctuations with STAR data [|3]. 
The curves and data are for beam energies 20, 130, and 200 GeV. The behavior of 
^dyn appears inverted because this quantity is negative. 



I find that Vdyn and (5ptiSpt2) I ' {Pt) 2 minimize trivial acceptance effects and are 
correspondingly less ambiguous. 
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